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Abstract 

We discuss, within the framework of the Standard Model, the calculation of the two- 
loop electroweak contributions to the anomalous magnetic moment of the muon involving 
triangle fermionic loops of leptons and quarks. Because of the large ratios of masses 
involved, these contributions are rather large. The result we obtain differs from a previous 
estimate reported in the literature. The discrepancy originates in the cancellation of 
anomalies in SU(3) C x SU(2)l x U(1)y, a cancellation that requires the consideration of 
both leptons and quarks within each generation and that had been previously overlooked. 
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1. There is a forthcoming experiment at the Brookhaven National Laboratory which plans 
to measure the anomalous magnetic moment of the muon with an expected uncertainty of 
±40 x 10~ n . This will correspond to an improvement by a factor of twenty with respect to 
the latest result obtained from the experiment performed at CERN, which gave (T|: 

a» = jfoi " 2 ) = 11659230(85) x 10~ 10 . (1) 

The muon g — 2 project at BNL M has renewed the interest on the part of some theorists 
to improve the accuracy of the corresponding prediction in the Standard Model. Here we 
want to concentrate on the electroweak higher-order corrections. The history and the present 
status of the calculation of the dominant electromagnetic contributions to can be found 
in the series of review articles ||, [Q], and 1J. A recent phenomenological re-evaluation of 
the hadronic vacuum polarization effect on has been made in ref. ||. The hadronic light- 
by-light scattering effect on has also been recently reconsidered, within the framework of 
low-energy QCD, in refs. {|, || and §§. 

The one-loop contributions to due to the electroweak interactions of the Standard 
Model were calculated quite a long time ago []!(]] . The relevant Feynman diagrams are shown 
in Fig.l. The corresponding result is 



w , ak _ Gp f 10 4 2 _ 2 / m£ Mf \ yi x\2 - x) 1 



where and a M are the vector and axial-vector couplings of the Z to the muon. In general, 
for a fermion / , 

v f =lf -2Q f sm 2 e w , a f = lf. (3) 

The contribution from the Higgs, given in terms of a parametric integral in eq. decouples 
in the infinite mass limit. Numerically, with the Higgs contribution neglected, 

a Weak = 1Q5 x 10 -11 _ (4) 

The leading two-loop electroweak contributions to have been discussed in ref. 
These authors have selected all the possible sources of logarithmically enhanced terms which 
appear because of the large ratios of masses involved, and which result in contributions to 
of order 

with M 3> m, where typically M is the Z mass and m a fermion mass, like for instance the 
muon. Apart from these, other possible logarithms involving the Higgs mass like log MnjMz 
or log Mff/mt might also appear but, unless Mh 3> Tn t , they are not so large and are dis- 
regarded. The authors of ref. [11] thus find the following overall correction Ao^ eak to the 
one-loop result in eq. (||): 

Aa^ eak ~ -42 x 10" 11 ; (6) 
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i.e. a rather large negative correction, of the same size as the expected experimental uncer- 
tainty hj. 

2. We shall be concerned with a specific class of the two-loop electroweak contributions: 
those induced by virtual fermionic triangle loops, represented by the Feynman diagrams in 
Fig. 2. The authors of ref. (T^j have only considered the subclass of these contributions where 
the fermion in the triangle loop is a lepton. In the 't Hooft-Feynman gauge and keeping only 
the asymptotic contributions from the large ratios of masses involved, the results they find 
are: 
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-9.11 x 10 
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-4.77 x 10 
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(7) 
(8) 
(9) 



One should realize that the results of eqs. (0-^) are actually gauge dependent. These 
results, as they stand, stem from the g^ u part of the Z propagator in the diagrams of Fig. 
2. 0. In the unitary gauge, for instance, the k^k u piece in the Z propagator yields an extra 
contribution to eqs. (0-^) that has a part that is common to these three equations. This extra 
common contribution is actually divergent and originates in the anomaly that results when 
one multiplies the triangle by k^. It is only when one sums over a complete generation that 
the anomaly vanishes and the result is finite and gauge invariant.^ As a consequence, strictly 
speaking, only the contribution of a full generation may be considered physically meaningful, 
but not that of a single fermion. 

The Feynman diagrams in Fig. 2 where the fermion in the triangle loop is a quark cor- 
respond to a new class of hadronic contributions to the muon anomaly. We shall call them 
the hadronic Z-7-7 contributions. The fact that the fermionic triangle subdiagram in Fig. 
2 has an Adler, Bell-Jackiw VVA anomaly, which in the Standard Model cancels when all 
the fermions of the same generation are included, implies the vanishing of the whole triangle 
diagram in the limit of exact mass degeneracy within each generation. We therefore expect 
important cancellations within each generation, which questions the overall estimate in eq. 
(|TJ) quoted from ref. [p|[ . 

We propose to investigate this problem first within the framework of effective quantum 
field theories where the underlying physics can be easily understood. We have also done an 
exact calculation of the contribution to from the diagrams in Fig. 2, and checked that, 
in the appropriate cases, the various asymptotic limits reproduce the simple effective field 
theory expressions. We have reproduced in particular the results in eqs. (0), ©, and (|9|) in 
the 't Hooft-Feynman gauge. We shall then discuss the new numerical results. 



3. Let us first consider the limit where the Z-mass is much larger than the masses of any 
of the other particles involved, which is certainly the case for the first and second generations, 

*The associated diagrams with a would-be Nambu-Goldstone have an extra suppression due to the mass of 
the fermion going around the triangle, m/, that goes like mi/M%. 

^The gauge invariance of the Standard Model is a fact only after the cancellation of anomalies is effected. 
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and keep only the leading log Mz contributions. When the Z-field in the Standard Model 
is integrated out, there appear new local four-fermion couplings induced by the tree-level 
exchange of the underlying Z propagator. These four-fermion couplings lead to two-loop 
diagrams like the ones shown in Fig. 3, which yield contributions to the muon g — 2 that 
are logarithmically divergent. This divergence is to be interpreted as cut off by Mz- If 
furthermore the fermion masses in the triangle loop are neglected with respect to the muon 
mass, the result from each fermionic contribution is the same, up to a factor Q 2 jdf proportional 
to the square of the electric charge Qf of the fermion times its axial coupling a t defined in 
eq. For quarks, there is an extra factor of three from the number of colours. The overall 
result we find for the muon g — 2 from the first generation, in the limit where m e = 0, and in 
the chiral limit where m u = rrid = 0, is then: 

A n „^ a -EQh, + £ + < 10 > 

where the 1 /e term encodes the ultraviolet divergences in dimensional regularization and the 
my-independent constants A, B are gauge dependent. They vanish in the 't Hooft-Feynman 
gauge but not in the unitary gauge. Each fermion contribution is separately gauge dependent 
and it is only the sum over the full generation that is physically meaningful. The fact that 
£/ Q)a f = is of course a property due to the anomaly cancellation in the SU(3) C x 
SU(2)l x U(l)y gauge theory. We shall later come back to a more elaborate discussion of the 
contribution from the first generation, which takes into account the effect of hadronic mass 
scales in the light quark sector of QCD. 

4. The evaluation of the leading contribution from the second generation of fermions is a 
little more delicate. We can still consider the effective field theory where the Z-field has been 
integrated out, but now the fermion masses in the triangle loop cannot be neglected with 
respect to the external muon mass. The contribution from the strange quark, in particular, 
requires a special discussion depending on whether or not one is willing to consider the 
chiral SU(3) limit where m s = 0, and on whether or not one wants to discuss spontaneous 
chiral symmetry breaking effects and the effects due to the QCD U(1)a anomaly. To a first 
approximation, we shall also consider the chiral limit for the strange quark and, as in the 
case of the first-generation estimate, we shall neglect for the time being the effect of hadronic 
mass scales. We then find : 



, w Pa w, GpmUf, L Ml 1, I , Ml 4,1 , Ml) 

= *^|«41„g 4.-5.4X10-", (11) 
a/2 8tt z it 

where the numerical result is for m c = 1.3 GeV. As expected on first principles, the Z-mass 
does not appear in the final result. 

5. The evaluation of the leading contribution from the third generation brings in an 
interesting issue, related to the fact that the top quark is heavier than the Z. Within the 
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framework of effective field theories, we now have to consider the case where the top field is 
integrated out first, corresponding to the limit mt 3> Mz- In this limit, the top quark in the 
Z-7-7 vertex decouples. The corresponding effective local Z-7-7 coupling induced by the 
top triangle loop goes as 1/ m 2 and induces a contribution to the muon g — 2 at the one-loop 
level via the Feynman diagrams shown in Fig. 4. If we work in the limit m t —* 00, this 
contribution vanishes. In this limit and in the unitary gauge (where there are no would-be 
Nambu-Goldstone bosons) the integration of the top leaves no trace behind in the form of 
new local effective operators relevant to the gyj — 2.^] 

In the effective Lagrangian without top the part of the Z propagator yields for r and 
b in the loop a finite contribution completely analogous to that found in the previous case 
for the first two generations. The k^k u part, on the other hand, yields the anomaly when 
multiplied by the triangle and this contributes a logarithmically divergent quantity that is to 
be interpreted as cut off by the top mass. Gathering all the pieces we obtain: 

e =T,b 6 

(12) 

where we have used m T = 1.78 GeV, m& = 4.3 GeV, and mt = 170 GeV for the numerical 
estimate. The last term in this equation comes from the k^k v part of the Z propagator. This 
last term with the sum over r and b would have been the contribution from the top quark 
had we done the calculation with the full theory, i.e. without integrating out the top. 

Altogether, treating the light u, d, and s quarks in the chiral limit, and neglecting the 
effect of hadronic mass scales, we find that in the Standard Model the leading contribution to 
the muon g — 2 from the full set of fermionic triangle graphs in Fig. 2 represents a correction 
of ~ —7% to the dominant one-loop electroweak contribution. 

Of course the limit mt —* 00 is quite far away from the real situation, and in the next 
section we will present exact expressions that will give corrections to this extreme limit. 

6. The effective field theory framework is not very useful when non-leading effects in the 
large ratios of masses are also required. In that case it is simpler to do a direct calculation. 
The result corresponding to the diagrams in Fig. 2, in the unitary gauge, for a fixed fermion 
/ in the triangle loop and without approximations has the form 

,„ , G-w ml a n ( m2 f M% m 2 f m 2 f M% 1 

AaJ^\ f = -Q 2 f a f - , (13) 

M y/28ir 2 it 1 1 { l m 2 m 2 Af§ l m 2 my) 

where we have separated the contributions generated by the g^ u term in the Z propagator (the 
^-"-function) from those generated by the k^k v term (the (/-function). As already mentioned 
earlier, it is the k^k u term which also generates a divergent piece, independent of the fermion 
mass in the loop, and which here has been subtracted. (It cancels of course, once the sum over 
the fermions of a generation is made.) The functions J- and Q have rather compact Feynman 
parametric representations. With x and y the Feynman parameters of the triangle loop, and 
u,v, and w those of the other loop we have: 

^There are of course other operators relevant to other processes. 



1 G-w m,, a ( Ml M% 
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-9 x KT 11 
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and 



r m 2 f M% /"I /"I /"I /"I fl-a! 8U 2 V 
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The contribution from the (/-function to the muon g — 2 is only sizeable for the t quark. 

2 

It is in fact this piece which reproduces the leading log j^r discussed earlier. We can now 
obtain the non-leading behaviour as well, with the result: 



a Weak mi _ G F m\ a (8 m 2 16 2Mf a /Mf\l 



(16) 



We can also use the expressions above to obtain a more accurate evaluation of the con- 
tributions to Aa^ eak when the fermion / in the triangle loop is one of the three leptons : 
/ = e, fi, t. When the fermion / in the triangle loop is a quark, the expression in eq. ( |l3| ) 
corresponds to the evaluation one obtains in the limit where the QCD gluonic interactions 
are neglected. If the quark in the triangle loop is a heavy quark: f = c,b,t, the QCD gluonic 
interactions can be treated perturbatively, and the corrections to the lowest-order estimate 
obtained, will be down by a typical factor a s (fi 2 )/ir, where the mass scale fi runs between m 



f 



and Mz [13|. We can therefore use eq. ( |i3|) to reliably obtain a more elaborate estimate of 



the contribution to the muon g — 2 from the third generation of fermions, with the result : 

AaJ^Mlkfls) = "8.2 x HT 11 , (17) 



to be compared with our approximate estimate in eq. fll2j). 

7. When the fermion / in the triangle loop is a light quark: f = u,d, s, the result obtained 



from a straightforward application of the expression in eq.(13) with either current algebra 
quark masses or constituent quark masses can be rather misleading. QCD perturbation 
theory is not justified in this case; and it would be erroneous to use it since, among other 
things, it neglects the fact that chiral symmetry is spontaneously broken. An appropriate way 
to discuss this problem is within the combined framework of chiral perturbation theory (xPT) 
and the 1/N C expansion. (For a recent review, where earlier references can also be found, see 
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e.g. ref. ||T4| .) To lowest order in the chiral expansion in U(3)l x U(3)r, the hadronic Z-7-7 
interaction appears via the one Goldstone meson exchanges between the effective coupling 



2 sim?w cos i?w V V3 V6 / 
induced by the lowest 0(p 2 ) chiral effective Lagrangian, and the effective 0(p 4 ) coupling 

i.e. the term in the Wess~Zumino Lagrangian which reproduces the Adler, Bell-Jackiw 
anomaly. The corresponding Feynman diagrams are shown in Fig. 5. The evaluation of the 
contribution to the muon g — 2 from these diagrams, in the unitary gauge and in the chiral 
limit, leads to the result : 



3 1Og 4 + 9 + O (4 1Og 4)J= 5 - 0Xl ° ' (20> 

where we have also subtracted the divergent piece generated by the k^k u term in the Z 
propagator, which cancels with the corresponding pieces generated by the e, /i, and c fermionic 
loops. We wish to emphasize that to lowest non-trivial order in xPT, in the chiral limit and in 
the large N c limit of QCD, this is an exact result. When added together with the contributions 
to Aa^ ca,k \f in eq. ( |l3| ) from the leptons e, [a, and the c quark in the triangle loop, it gives as 
the result to the contribution to Aa^ eak from the first and second generation of leptons and 
quarks : 

Aa)f eak [e,(f,u;/i 1 s,c] = -8.7 x KT 11 , (21) 
to be compared with our first estimate, the sum from eqs. (|To| ) and 



Aa^ k [u,d,s] 



J2 87T 2 7T 



XPT 



8. We now want to discuss the sources of corrections to the above calculation and their 
possible estimate. 

One source is the effect of the a s (/j, 2 )/TT perturbative corrections to the lowest-order c, 
b, and t loop calculations which we have already mentioned. They are, in principle, under 
control and one does not expect them to drastically change our results. 

The other sources of corrections concern the light-quark sector. Here, besides the obvious 
change in the Nambu- Goldstone propagator due to finite pseudoscalar masses, we expect chiral 
loop corrections due to the explicit breaking of the chiral symmetry, as well as corrections 
due to the contributions from higher order terms in the effective chiral Lagrangian. Chiral- 
loop corrections are suppressed in the 1/N C expansion. We expect them to give corrections 

C(jg^rp)) perhaps enhanced by chiral logarithmic factors. However, the most important 
effect to next to leading order in the 1/N C expansion is likely to be the fact that, because of 
the U{1)a anomaly, the singlet t)q particle acquires a large mass. The effect of this mass will 
be to damp the contribution of the rj to the muon g — 2 in eq. (|20|). 

Concerning the effect due to higher order terms in the chiral expansion, it is possible to 
make estimates using models of the QCD low-energy effective action at large -/V c , which have 
been developed during the last few years (see e.g. ref. The simplest version of these 
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models amounts in practice to giving a constituent mass Mq to the u, d, and s quarks, and 



to modulate the axial- vector coupling of the constituent quarks with a constant qa Hlgl , [16|. 
Here, care must be taken, however, on the way constituent quark masses and the coupling g& 
are introduced. As already mentioned, the VVA vertex in the triangle loop in Fig. 2 has an 
anomalous Ward identity : in the chiral limit the VV<9A vertex has a universal form which 
is, in particular, at the origin of the effective chiral realization discussed above. With k^ the 
momentum flowing from the axial coupling, and T^ a ^ the full VVA vertex, we can decompose 
this vertex as follows : 

T"°0 = ^k^kpTP^ + (t^ - ^k^kpTP^ = R^^j . (22) 

The first term in the r.h.s., when modulated by the appropriate X)/Q/ a / factor, with 
f = u,d,s and including the iV c colour factor reproduces, in the chiral limit, the expres- 
sion obtained in the effective chiral realization that we have already discussed (see eqs. ( |l8| ) 
and ©.) This is an example of the 't Hooft anomaly matching condition |l7| applied to 
QCD. The constituent chiral quark picture should then be applied to the second term only, 
i.e. the R^ a ^ tensor defined as the difference in eq. (p2|). (Notice that the contribution 
to the muon g — 2 from the .R-tensor is gauge invariant.) This way, and using the values 
Mq ~ (265 ± 10) MeV and gj± ~ 0.6 ± 0.1, we find a correction of at most ~ 30% to the result 
in eq. (|20|). From this we conclude that 50% is a safe estimate of the size of the expected 
errors in the xPT calculation in eq. (|20|). 

Our final estimate of the leading (in the sense of eq. (5)) electroweak contributions from 
the full set of fermionic triangle loops in the Standard Model, from eqs. (|l7j) and (|2l|) is then : 



Aa^ eak [e,(i,u; jj,,s,c; r,b,t] = -(16.9 ±2.5) x 10 -11 . (23) 
This result, when added to the other leading two-loop electroweak corrections calculated 



in [11 1, corresponds to an overall estimate 



Aa^ eak ~ -(36.9 ± 2.5) x 10~ n ; (24) 
i.e. a negative correction of ~ 19% to the lowest order electroweak contribution in eq. (Q). 

S.P. would like to thank E. Alvarez, L. Alvarez-Gaume and S. Yankielowicz for conversa- 
tions. We would like to thank A. Czarnecki for pointing out a sign error in the contribution 
of the top in a previous version of the manuscript. 



Note Added 

After completion of this work we have become aware of a paper by A. Czarnecki, B. Krause 
and W.J. Marciano, |hep-ph/9506256 , where a subset of the constant (i.e. non-log M/m 
enhanced, see eq. (5)) contributions has been calculated. These contributions originate 
mainly in m^-like corrections; they are negative and push the final result in eq. (24) to 
-45 x 10 -11 . 
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Figure Captions 



Fig.l : Feynman diagrams which, in the Standard Model, give the lowest order electroweak 
contribution to the muon g — 2. 

Fig. 2 : Two-loop electroweak corrections induced by virtual fermionic triangle loops. 

Fig.3 : Two-loop electroweak contributions, with virtual fermionic triangle loops, induced by 
a local four-fermion effective coupling. 

Fig.4 : One-loop electroweak contribution induced by an effective Z-7-7 local coupling. 

Fig.5 : Lowest order contribution in the effective chiral realization of QCD, which leads to a 
Z-7-7 hadronic contribution to the muon g — 2. 
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Fig. 1 




Fig. 3 
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Fig. 5 



